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Abstract 



At a certain finite neck radius during the descent of a fissioning nucleus from the 
saddle to the scission point, the attractive nuclear forces can no more withstand 
the repulsive Coulomb forces producing the neck rupture and the sudden absorp- 
tion of the neck stubs by the fragments. At that moment, the neutrons, although 
still characterized by their pre-scission wave functions, find themselves in the newly 
created potential of their interaction with the separated fragments. Their wave func- 
tions become wave packets with components in the continuum. The probability to 
populate such states gives evidently the emission probability of neutrons at scis- 
sion. In this way, we have studied scission neutrons for the fissioning nucleus 236 U, 
using two-dimensional realistic nuclear shapes. Both the emission probability and 
the distribution of the emission points relative to the fission fragments strongly de- 



pend on the quantum numbers of the pre-scission state from which the neutron is 
emitted. In particular it was found that states with = 1/2+ dominate the emis- 
sion. Depending on the assumed pre- and post-scission configurations and on the 
emission-barrier height, 30 to 50% of the total scission neutrons are emitted from 
1/2+ states. Their emission points are concentrated in the region between the newly 
separated fragments. The upper limit for the total number of neutrons per scission 
event is predicted to lie between 0.16 and 1.73 (depending on the computational 
assumptions). 
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1 Introduction 



When a cold nucleus undergoes fission, the coupling between the collective and 
intrinsic motions can either be adiabatic (the nucleons remain in the lowest 
energy levels) or diabatic (the nucleons are promoted to excited states). This is 
mainly governed by the velocity of the collective motion relative to the Fermi 
velocity. 

At the beginning of the descent from saddle to scission, the Coulomb repulsion 
is compensated to a large extent by the nuclear attraction. The coupling be- 
tween the collective and intrinsic motions is therefore weak (quasi-adiabatic): 
the collective motion is relatively slow and only few nucleons are excited with 
increasing deformation. 

Close to scission however this quasi-equilibrium breaks down. The neck be- 
tween the fragments ruptures and it is quickly absorbed by the fragments. The 
coupling is then diabatic. 

In the present work, we are making the extreme assumption that the fission 
motion is adiabatic until the neck ruptures at finite radius and that this rup- 
ture is followed by the sudden absorption of the neck protuberances. At that 
moment a transition between two quite different nuclear configurations occurs. 
The corresponding large difference in the potential energy of deformation is 
made available to the internal degrees of freedom, that is, mainly to the nu- 
cleons but also to preformed clusters. 

This sudden release was proposed as an emission mechanism by Halpern [1] 
without being quantitatively developed. The only attempt so far to carry out 
calculations using the sudden approximation was done later using a simplistic 
one-dimensional model [2J for the emission of scission alpha particles. A grad- 
ual potential change and its effect on the energy transferred to the nucleons 
was studied by Fuller [3] and by Boneh and Fraenkel [I]. 
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In the present paper we have generalized the formalism from Ref. [2] to two 
dimensions in order to include realistic nuclear shapes, and have applied it to 
the emission of neutrons during 235 U (n^,f) reaction. 

Emission of scission nucleons has always been an intriguing subject. Based 
on energetic considerations they should be the most probable among the light 
particles that accompany the fission process. In particular they should outnum- 
ber the alpha particles. Experimentally however the emission rate for scission 
protons represents only few percent of that for scission alphas [5j while the ex- 
isting results on scission neutrons are very contradictory 05]. Understanding 
this paradox is a central piece of the nuclear fission puzzle. 

For scission neutrons, there is a clear need for more precise experimental data. 
Nevertheless, in a recent paper, Kornilov [7] has re-analyzed three indepen- 
dent experiments on 252 Cf(sf) with a more realistic model for the description 
of the neutron evaporation from the fragments. In all these experiments, the 
neutron energy and angular distributions relative to the fragment direction 
were measured. Kornilov showed that a good agreement exist between these 
experiments, and that a neutron surplus of (30 ±5)% exists at about 90° 
relative to the direction of the moving fragments. These neutrons do not orig- 
inate from fully accelerated fragments, and would represent about 10% of the 
total fission neutrons yield. A similar conclusion has been drawn from angular 
distribution of neutrons measured by Franklyn [8] for the 23S U(nt^,f) reaction. 
This angular distribution could not be reproduced if scission neutrons are ne- 
glected but could be very well fitted with a scission neutron contribution to 
the total neutron yield of 20%. 

It is also worth noting that using a systematics of ternary charged particle 
yields established by Val'ski [9], scission neutron yields can be extrapolated for 
various fissionning nuclei. These scission neutron yields are found to be (0.18 ± 
0.04)/fission for 252 Cf(sf) and (0.55 ± 0.09)/fission for 235 XJ(n th ,f), representing 
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about 4.8% and 22.7% of the total fission neutrons yields, respectively. A 
recent review of experimental data on the scission neutron characteristics can 
be found in [TU] . 

In the next section, we will show how to calculate primary fragments excita- 
tion energies, scission neutrons multiplicities and spatial distributions of the 
emission points starting from the single-neutron wave functions at the two nu- 
clear configurations between which the sudden emission is supposed to occur. 
Section 3 contains our predictions for the low energy fission of U-236 obtained 
for independent and pairing correlated neutrons, for two sets of initial and 
final configurations, and for two energy cuts separating bound from unbound 
states. In section 4, we will give a brief summary of the main results. 

2 The formalism of the sudden approximation 

The single-particle wave functions for an axially-symmetric fissioning nucleus 
have the general form 

(p, z, 4>) = u(p, z)e< a ~*)*\ T> +d{p, z)e< n+1 i>\ |>, (1) 

where u(p, z) and d(p, z) contain the spatial dependence of the two compo- 
nents, spin up and down respectively. Q is the projection of the total angular 
momentum along the symmetry axis and is a good quantum number. Since we 
are dealing with symmetric fission, the parity 7r is also a constant of motion. 

If the scission is characterized by a sudden change of nuclear deformation, an 
eigenstate of the "just before scission" hamiltonian will be distributed over 
the eigenstates of the "immediately after scission" hamiltonian: 

\¥>= £ a if \^>, (2) 

allstates 
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where 



--< ^ f \¥ >=2nj J (uV + d l d f ) pdpdz. (3) 



One can notice that only |^ > states with the same (0, 7r) values as > 
will have non-zero contributions. These states are mainly bound states but 
contain also a few discrete states in the continuum that will spontaneously 
decay. Therefore the emission probability of a neutron that had occupied the 
state 1^ > is 

pL= E k/l 2 = i- E k/l 2 - ( 4 ) 

unbound-states boundstates 

In this way we do not need to use the states in the continuum that are less 
precise in the numerical diagonalization. The Pl^s will be referred as partial 
probabilities since they only consider one occupied state. The limit between 
bound and unbound states is the barrier for neutron emission, which is zero 
if the centrifugal potential is null. At this point it is worth noticing that our 
criterion for emission is based on energetic considerations, while Halpern's was 
based on spatial considerations. 

Summing over all occupied states one obtains the total number of scission 
neutrons per fission event 

Nn = T, V l P L, (5) 

i 

where vf is the ground-state occupation probability of >. 

We have also calculated the part of the initial wave function that was emitted 



\¥ em >=\¥>- £ a if \^>, (6) 

boundstates 
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since it gives access to the spatial distribution of the emission points 



Sem{p, Z) = \% m {p, z) 



(7) 



Of course, we have 



J S em (p, z)pdpdz = N n . 



(8) 



Finally, the occupation probabilities after the sudden transition are obtained 
by summing the contributions of all initial states that populate a given final 
state: 

v? = E«?l«d 2 - (9) 

i 

They are of interest since they show the degree of excitation in which the 
fragments are left: 

E FF = E V}e*- £ v}ef, (10) 

boundst. boundst. 

where e* are the energies of the bound eigentstates of the hamiltonian "im- 
mediately after scission" . 



3 Numerical results 



We have applied the formalism described above to quantitatively study the 
emission of scission neutrons during the low-energy fission of 236 U. 

The nuclear shapes involved were chosen to be Cassini ovals since they require 
only one deformation parameter e to describe realistic scission configurations. 
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To see how critical is the choice of the two deformations between which the 
sudden transition is supposed to occur, we have used two sets of values (e*, e/), 
0.975 — > 1.010 and 0.985 — > 1.001, describing a longer and a shorter jump 
respectively (e—1 describes a zero-neck scission shape). An idea of the shapes 
involved is given (Fig. [1]) by the neutron equipotential lines corresponding to 
Vo/2, Vo=-40.2 MeV being the depth of the Woods-Saxon interaction. The 
neck radii at rupture are 1.9 fm and 1.5 fm respectively. 

The corresponding hamiltonians have been diagonalized in a deformed os- 
cillator basis and complete sets of eigenfunctions |\P*(ej) >, (e/) > and 
eigenvalues e l (ej), e* (e/) have been obtained [11]. From these, only the bound 
states (e < 0) enter in the present calculations. It turned out that 236 U has 98 
such states at all four considered deformations. 

3.1 Partial emission probabilities 

For each initial bound state we have calculated the expansion coefficients 
dif, Eq. (j3J), corresponding to all final bound states. Using these coefficients, 
we have calculated the partial emission probabilities, Eq. (j4j). The results are 
presented in Figures [2] and [3] for the transitions 0.985— »1. 001 and 0.975— >1. 010 
respectively. They are separated according to their quantum numbers. 

As expected the longer jump produces more scission neutrons. One can also 
notice that most of these neutrons are emitted from Q = 1/2 states with 
a preference for 7r = +• The variation of P l em with the energy of the single- 
neutron state is different for deep lying states than for states around the Fermi 
energy (ep ~-5 MeV). At the beginning (up to about -8 MeV), there is a slow 
increase; afterwards, the increase becomes steeper and attains a maximum 
(at about -4 MeV) followed by a moderate decrease. Although this type of 
variation is present for all fi-values, the absolute values of the partial emission 
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probabilities are decreasing with increasing fi. 

Over this average trend there are fluctuations due to the structure of the 
neutron states. To understand these features one has to take a closer look at 
the neutron wave functions involved. 

Figure H] shows two pairs of deep initial states together with their distribution 
of probabilities over final bound states. A pair is defined by two states that 
only differ by parity. The sudden transition was supposed to occur between 
= 0.985 and £/ = 1.001. A very low partial emission probability is asso- 
ciated with an initial wave function that is essentially distributed only over 
bound final wave functions. In general, this situation occurs when the initial 
wave-function resembles very much its adiabatic-partner final wave function. 
The wave packet in the final configuration is very narrow and the neutron re- 
mains essentially bound. Negative-parity initial states can more easily satisfy 
this condition since they cancel in the origin as do the states of the newly sep- 
arated fragments. Hence they are already in a post-scission configuration, and 
this explains why they are less emitted than the corresponding positive-parity 
states. For a better understanding of this last remark, we have compared in 
Fig. [5] a pair of wave functions at e;=0.975 with the corresponding adiabatic 
pair at e/=1.010. 

Figure M shows a case in which the nuclear structure effects are strong enough 
to invert the average trend. The states that are more bound (-15.3110 MeV 
and -15.9955 MeV) have a higher emission probability than the states that 
are less bound (-12.2621 MeV and -12.4234 MeV) since, as it can be seen, the 
latter wave packets are narrower. 

Figures [2] and [3] show that Pl m {Q = 1/2) exhibits a maximum at about - 
3 MeV. When the energy of the initial state increases there are fewer bound 
states with higher energy left in the final configuration and the distribution 
of |a;/| 2 becomes asymmetric, i.e. the higher-energy part of the distribution 
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is cut at zero (see right-hand side of Fig. [7]). Therefore the wave packet has 
to include more unbound states and this leads to an increase in the emission 
probability. However, loosely bound states (-0.3173 MeV and -0.7885 MeV) 
lie too high to be much perturbed by the sudden increase of the potential in 
the neck region (see Fig. [7]). It can be seen that, similar to the strongly bound 
states (Fig. 0J these states are essentially distributed over only one final state 
and have therefore a smaller partial emission probability. Moreover such states 
are practically unoccupied (vf ~ 0). 

Figure [H shows why states with higher Q- values (5/2) are less emitted. The 
presence in the deformed potential of a term in (Q ± 1/2) 2 / p 2 (reminiscent of 
the centrifugal potential) pushes the neutron away from the z-axis. Therefore 
their presence in the neck is reduced even for it = + states. Of course states 
with the highest projections of the angular momentum (Q = 9/2 and Q = 
11/2) lie very far from the z-axis and they can only be located in the future 
fragments (see Fig. [9]). Consequently, their emission probability is negligible. 

To see the effect of a larger difference between the "just before" and "immedi- 
ately after" scission configurations, we have plotted in Fig. [10] the same states 
as in Fig. @]for the transition 0.975 — > 1.010. As expected, the values P\ m are 
larger (by a factor from 4 to 10) since the coefficients |aj/| 2 are now distributed 
over more final states. 

3.2 Scission neutron multiplicities 

Summing the partial neutron emission probabilities P % em of all initially bound- 
states weighted by their occupation probabilities vf, Eq. El we have obtained 
the total scission-neutron multiplicities N n . The values obtained with several 
computational assumptions are found in Tables [1] and [2j 

For the dependence of vf on the particular state i, we have considered the 
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neutrons 



(1) independent: 




1 for i< 144/2, 



(11) 



otherwise 



(2) pairing-correlated: 




1 - 



e; - A 



(12) 



^(e, - A) 2 + A 2 



The pairing gap A and the Fermi energy A are the solutions of the two 
coupled-nonlinear BCS equations [12J. The closely related smooth gap pa- 



mass differences in nuclei throughout the periodic table. 

Tables Q] and [2] correspond to different separation lines between bound and 
unbound states in Eq. HI and 2 MeV respectively. The first value means 
that the centrifugal barrier is null. This is adequate when considering Q = 1/2 
states, since their /-distribution is highest at / = 0. However, states with higher 
f2-values are influenced by the 2D centrifugal barrier expressed by the term 
(f2±l/2) 2 /p 2 present in the hamiltonian. Eigenstates with energies lower 
than the top of the centrifugal barrier are quasi-bound states, and could there- 
fore be considered as still present in the fragments, their reabsorption by the 
fragments being large. The 2 MeV results should be viewed as an indication of 
the sensitivity of the results on this cut-off energy parameter. This sensitivity 
is rather important since the total emission probabilities are roughly reduced 
by half. 

Due to the uneven distribution of states with given Q, the variation of N n 
with the energy cut-off is not monotonic. This explains why in the last three 



rameter A ~ 13.3/ y/A MeV was taken from the systematics [13J of odd-even 
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columns of Table HJ N n is higher for 5/2 than for 3/2 states, contrary to the 
general behaviour of N n which is to decrease with Q. 

Analyzing these tables, one can also notice that more than 93% of the scission 
neutrons are emitted from 1/2, 3/2 or 5/2 states, more than 55% are emitted 
from 1/2 states and more than 33% from 1/2+ states. Although for most Q 
values, N n is higher for positive than for negative parity states, Q=3/2 is an 
exception. For those states the wave functions with negative parity appear to 
have a maximum in the neck region. This is possible since what counts for the 
behaviour is the parity of the functions d(p, z) and u(p, z) and not of the total 
wave function \I/ (see Eq. [TJ. For instance, 



is negative if u(p, z) is symmetric with respect to the p-axis. It turns out that 
u is the dominant component, i.e., J u 2 > J d 2 . 

Tables [T] and [2] contain also a comparison between the two sets of (e, — > ef) 
values considered. As expected, the calculated scission-neutron multiplicity 
increases with the difference (e/ — e^). We found a factor of about 5 for 0-MeV 
cut-off and of about 6 for 2-MeV cut-off. The sensitivity of our results on the 
choice of 6j and ej values is therefore very high, while the difference in neck 
radii between the two configurations is not so large (0.4 fm). 

3.3 Distribution of the emission points 

To obtain a detailed picture of scission neutron emission and to follow their 
motion thereafter one needs to know the distribution of the emission points. 

For this purpose, we have plotted in Figs. [TTJ C2] and [13] the square moduli 
of the emitted part of the neutron wave functions (Eq. E]) for the same states 




(13) 
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as in Figs. HI [9] and [TIB Since they represent the high-energy tails of the wave 
packets in the final potential well, they resemble the corresponding initial wave 
functions but they have more nodes (more structures in the (p, z)-plane). This 
is particularly clear for the highest Q values: 9/2 and 11/2 (Fig. [12]) . Because 
of the low level density of these states, the emitted part coincides with the 
next state in the spectrum, i.e., it has one more node than the initial state. 

Those wave functions that are initially present in the neck region will have 
higher emission probabilities. As predicted by Halpern [TJ, the wave function 
of the emitted neutron mainly contains the part of the initial wave function 
that was located in the neck. 

Summing over all initially-bound states like in Eq. (j7|), we have obtained the 
total distribution of the emission points (Fig. (T3j) for the two sets of 'just be- 
fore' and 'immediately after' scission configurations. Apart from the absolute 
values, the two distributions are similar. This can be better seen if we look only 
along the p and z axes as in Fig. [T5j The large majority of the scission neu- 
trons are emitted in the region between the fragments. To escape the nuclear 
re-absorption, they have to move perpendicular to the fission axis. With much 
lower probability (three orders of magnitude) unbound neutrons are however 
present also inside the fragments and at the poles of the fissioning system. 
This latter component is either re-absorbed or emitted in the direction of the 
fission axis (polar emission). 

3-4 Primary fission fragments excitation energies 

A sudden transition at scission not only produces neutrons but also leaves the 
primary fragments in an excited state. This amount of excitation energy is 
important to know since it can be used to emit neutrons during and after the 
acceleration of the fragments. The later are called 'prompt' neutrons and have 
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been extensively studied both theoretically and experimentally (for a review, 
see [H] and references therein). 

Figure [16] shows the probabilities to occupy bound-states in the final configu- 
ration according to Eq. (Q. Comparing with the same probabilities in the BCS 
ground-state, one can see how states above the Fermi level have been popu- 
lated at the expense of states below. As for the emission probabilities (Figs. [2] 
and Fig. [3]), the occupation probabilities strongly depend on the structure of 
the individual states involved. These structure effects are more pronounced in 
the lower part which corresponds to the larger difference of the deformations 
between which the sudden transition occurs (0.975— >1. 010). 

Using these partial occupation probabilities Vj in Eq. (fTUj) . we have calculated 
the total excitation energy of the fragments immediately after scission. We 
found 12.2 MeV for the shorter jump and 51.88 MeV for the longer jump, i.e., 
roughly the same factor 4-5 as for multiplicities. 

4 Summary and conclusions 

We have presented a model calculation of quantities that characterize the nu- 
clear configuration at scission: excitation energy of primary fission fragments, 
scission neutron multiplicity and initial distribution of their emission points. 
The approach is relatively simple and the results are easy to interpret micro- 
scopically in terms of the quantum numbers of the states involved. 

The results turned out to be highly sensitive to the value of the minimum 
neck radius assumed (when the neck ruptures). A precise measurement of the 
scission neutron multiplicity could therefore determine this important dynam- 
ical parameter. On the other hand, if the minimum neck radius is 1.5 fm (a 
reasonable value) 15% of the total neutrons emitted during the symmetric 



13 



fission of U-236 are scission neutrons. They are therefore non negligible and 
they should be explicitly taken into account in any fission data analysis. 

Our model also provided an estimate of the excitation energy of the primary 
fragments (12 MeV), a deciding factor for the emission of prompt neutrons, 
unavailable from other sources. 

Due to the sudden approximation used, the values obtained represent upper 
limits and this fact has to be remembered when comparing with experimental 
data. 
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Table 1 

Scission neutron multiplicities with "0-MeV" cut-off. 





N n (0.975 -> 1.010) 


N n (0.985 -> 1.001) 




BCS 


fit.pn 'Fiinrtinn 


BCS 


St.PD 'Fiinrtinn 


1/2 -/+ 


n 371 /n ^fi7 


n 4^n /n ^8 


n nfi9 /n 144 


n nQ7 /n 1 4^ 


3/9 -/4- 
o/z /-\- 


D 914/fl D89 


u . loo u.uuy 


n D98 /n ni 3 


n m 7/n m n 


5/2 -/+ 


1 88/0 909 


930 /0 941 


049/0 043 


0^1 /0 0^9 


7/2 -/+ 


0.047/0.047 


n n ^7 / n n n ^7 

0.007/0.007 


0.011/0.011 


r\ r\r\i I r\ aai 

0.001/0.001 


9/2 -/ + 


0.006/0.006 


0.002/0.002 


0.001/0.001 


0.000/0.000 


11/2 -/+ 


0.001/0.001 


0.000/0.000 


0.000/0.000 


0.000/0.000 


Total -/+ 


0.828/0.904 


0.845/0.878 


0.144/0.213 


0.168/0.209 


Total 


1.733 


1.723 


0.358 


0.377 


Table 2 










Scission neutron multiplicities with "2-MeV" 


cut-off. 






N n (0.975 -> 1.010) 


iV n (0.985 -> 1.001) 




BCS 


Step Function 


BCS 


Step Function 


1/2 -/+ 


0.259/0.400 


0.306/0.381 


0.032/0.089 


0.045/0.085 


3/2 -/+ 


0.187/0.072 


0.139/0.062 


0.021/0.012 


0.014/0.009 


5/2 -/+ 


0.027/0.032 


0.030/0.033 


0.004/0.004 


0.005/0.005 


7/2 -/+ 


0.012/0.012 


0.002/0.002 


0.003/0.003 


0.001/0.001 


9/2 -/+ 


0.006/0.006 


0.002/0.002 


0.001/0.001 


0.000/0.000 


11/2 -/+ 


0.001/0.001 


0.000/0.000 


0.000/0.000 


0.000/0.000 


Total -/+ 


0.493/0.523 


0.480/0.481 


0.062/0.110 


0.065/0.099 


Total 


1.016 


0.960 


0.172 


0.164 
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Fig. 1. Equipotentials Vq/2 for the two scission configurations studied in this work. 
The radius of the neck at scission is 1.5 fm for the configuration (0.985,1.001) and 
1.9 fm in the case of (0.975,1.010). 
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Fig. 2. Partial emission probabilities P\ m for the configuration (ej, ej)=(0. 985, 1.001). 
The solid line represents the single-particle occupation probabilities in the initial 
configuration, calculated in the BCS formalism (right scale). The dashed line rep- 
resents the location of the Fermi energy at -4.926 MeV. 
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Fig. 3. Same as in Fig. [2] for the configuration (ej, e/)=(0. 975, 1.010). The value of 
the Fermi energy is now -5.229 MeV. 
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Fig. 4. Two pairs of low-lying initial states (left) and their expansion coefficients on 
the final bound states (right) for the configuration (0.985,1.001). Initial states with 
a high probability of presence in the neck are strongly emitted (large P l em ) since 
they can only expand on final states in the continuum. 
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Fig. 5. A pair of wave functions at = 0.975 with the corresponding adiabatic 
pair at e/ = 1.010. Initial states that are localized in the neck (upper graph) lead 
to large emission probabilities. On the contrary, states mostly present outside the 
neck (lower graph) are already in a post-scission configuration, and are much less 
emitted. 
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Fig. 6. Same as Fig. [4] for two different pairs of initial states. Although the states 
represented in the upper graph lie lower in energy than the ones in the lower graph, 
they have a higher emission probability, due to nuclear structure effects. 
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Fig. 7. Evolution of the emission probability spectrum as a function of the energy 
of the initial states. 
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Fig. 8. States with Q = 5/2 are less emitted than Q, = 1/2 states. The presence of 
the term ($7 ± l/2) 2 /p 2 in the deformed potential pushes the wave functions away 
from the z-axis and strongly hinders \ip\ 2 in the neck region. 
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Fig. 9. States with even higher fi-values (here, 9/2 and 11/2) cannot be present in 
the neck region. Therefore, their emission probabilities are negligible. 
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Fig. 10. The states corresponding to those of Fig. H] are represented here for the con- 
figuration (0.975, 1.010). As expected, the larger difference between the "just before 
scission" and "immediately after scission" configurations lead to higher emission 
probabilities. 
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Fig. 11. Emission distribution points corresponding to the initial wave functions 
represented in Fig. SI The solid line represents the equipotential Vq/2 for ey=1.001. 
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Fig. 12. Emission distribution points corresponding to the initial wave functions 
represented in Fig. [9j The emitted part of a state with high-f2 value corresponds 
almost perfectly to the same-0 value state lying just above in energy. The number 
of nodes is simply incremented by one unit. 
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Fig. 13. Emission distribution points corresponding to the initial wave functions 
represented in Fig. 1 101 for the configuration (0.975, 1.010). 
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Fig. 15. Variations of S em (p, z) along the z- and p-axes, for both scission configura- 
tions studied. 
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Fig. 16. Occupation probabilities of single- neutron levels "immediately-after-scis- 
sion" for the cases (0.985,1.001) (top) and (0.975,1.010) (bottom). 
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